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ABSTRACT 

We report on the redshift measurement and CO line excitation of HERMES J105751. 1+573027 
(HLSW-01), a strongly lensed submillimeter galaxy discovered in Herschelj 'SPIRE observations as 
part of the Herschel Multi-tiered Extragalactic Survey (HerMES). HLSW-01 is an ultra-luminous 
galaxy with an intrinsic far- infrared luminosity of Lfir = 1.4 x 10 13 L@, and is lensed by a massive 
group of galaxies into at least four images with a total magnification of \i = 10.9 ± 0.7. With 
the 100 GHz instantaneous bandwidth of the Z-Spec instrument on the Caltech Submillimeter 
Observatory, we robustly identify a redshift of z = 2.958 ± 0.007 for this source, using the 
simultaneous detection of four CO emission lines (J = 7 — > 6, J = 8 — >• 7, J = 9 — > 8, 
and J = 10 — > 9). Combining the measured line fluxes for these high- J transitions with the 
J = 1 — > 0, J = 3 — > 2 and J = 5 — >• 4 line fluxes measured with the Green Bank Telescope, the 
Combined Array for Research in Millimeter Astronomy, and the Plateau de Bure Interferometer, 
respectively, we model the physical properties of the molecular gas in this galaxy. We find that 
the full CO spectral line energy distribution is well described by warm, moderate-density gas 
with Tkin = 86 — 235 K and nn 2 = (1.1 — 3.5) x 10 3 cm~ 3 . However, it is possible that the 
highest-J transitions are tracing a small fraction of very dense gas in molecular cloud cores, 
and two-component models that include a warm/dense molecular gas phase with T^m ~ 200 K, 
riH 2 ~ 10 5 cm~ 3 are also consistent with these data. Higher signal-to-noise measurements of 
the J up > 7 transitions with high spectral resolution, combined with high spatial resolution CO 
maps, are needed to improve our understanding of the gas excitation, morphology, and dynamics 
of this interesting high-redshift galaxy. 

Subject headings: galaxies: high-redshift - galaxies: individual (HERMES J105751. 1+573027) - galaxies: 
starburst - submillimeter: galaxies 
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1. Introduction 

Galaxies selected at submillimeter (submm) 
and millimeter (mm) wavelengths (hereafter SMGs) 
are predominantly high-redshift (1.5 < z < 3.5) 
systems passing through an important, intense 
starburst phase of their evolution, with typical 
star format i on ra t es (SFRs) > 100 y r -1 (e.; 



Blain et all 120021: IChapman et all 120051) . Over 



the past 13 years, a large number of deep, wide- 
area surveys at A = 850 — 1200 (xm have been 



Scott et al.ll2002l iGreve et al. 2004: 


CoDDin et al. 


2006: iBertoldi et al. 


20071 Perera et al. 


20081 


WeiB et al.ll2009bl Austermann et al.ll2010l 


Scott et al. 


2010). These wavelengths sample the Ravleigh- 



Jeans tail of thermal dust emission, where this 
dust is heated by ultraviolet (UV) light from 
intense star formation or active galactic nuclei 
(AGNs). Owing to a strong negative fc-correction, 
galaxies with the same bolometric luminosity (a 
proxy for SFR) that are selected at A > 500 /im are 
equally detectable for 1 < z < 10 in a flux-limited 
survey. However, the relative ease of detecting 
a significant number of SMGs in deep, wide- 
area surveys is countered by the time-consuming 
multi-wavelength follow-up necessary to derive in- 
formation on their redshifts, SFRs and star for- 
mation efficiencies, morphologies, and dynamics, 
and on t he importance of AGNs in these sys- 
tems (e.g. 



Chapman et al.ll2005t IPope et al.| [200l 



Younger et al.ll2007l l2009t iBiggs et al.ll201C 

Measuring redshifts for a large number of 
SMGs is necessary for determining their con- 
tribution to the cosmic star formation h istory 
( Chapman et al. 2005 : Aretxaga et al. 20071 ). and 
is crucial for carrying out detailed studies of a 
representative sample of this population. Ob- 
taining redshifts through optical spectroscopy is 
difficult, given the poor positional accuracy of 
most SMGs (arising from the low resolution of 
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ground-based submm/mm telescopes) and their 
extreme dust obscuration. The most direct red- 
shift measurement for SMGs is through the detec- 
tion of the rotational transitions of carbon monox- 
ide (CO), since SMGs contain large reservoirs of 



molecular gas (~ IP 10 " 11 M Q ; iGreve et all 120051: 



Tacconi et al.ll2006l |2008( ). Until very recently - 
with the exception of three blind CO redshift de- 
tections (IWeifi et al.l l2009at ISwinba nk et alJl20ld 
Lestrade et alj[2010h - CO measurements of SMGs 
had been limited to sources with known, optically 



2004 



determined redshifts (e. g.lFraver et al.lll998l 199 
20081: llvison et al.ll200lt iDownes & Solomonll200_. 
Genzel et all 120031: iNeri et al.l l2003t ISheth et al 



2006 



Kneib et al.l2005l:lGreve et al.ll2005tlTacconi et ail 
Chapman et al.l 120081 iBothwell et al.ll2010t 



Engel et al.l l2010h . since the limited bandwidth 
and sensitivity of most receivers precluded a blind 
search for redshifted CO lines. 

Two recent developments have spurred the 
rapid growth of the number of SMGs with 
blind CO redshift detections over the past year. 
The first is the availability of sensitive, wide- 
bandwidth receivers designed specifically for the 
detection of one or more CO lines. These include 
Z-Spec, a grating spectrometer with 160 silicon- 
nitride micro-me sh bolometers ope r ating from 



190 - 310 GHz (Navlor et al. 2003: E arle et al 



avii 

20061 : iBradford et aTTT2009h . For z > 0.5, at 



least two CO transitions are redshifted into the 
Z-Spec bandpass, allowing for a robust redshift 
determination for galaxies at z < 3 (or higher, 
provided that the CO J up > 8 transitions are ex- 
cited). The second development is the recent de- 
tection of extremely bright SMGs with apparent 
far- infrared (FIR) luminosities > 10 14 L Q . Al- 
though a few discoveries of such extreme SMGs 
have been made in smal l -scale mapping sur - 
veys ( Swinbank et al.l l2010l : llkarashi et al.l 120101 ) , 
large area (10 — 200 deg 2 ) surveys carried out 
with the South Po le Telescope (SPT) at 1.4 mm 
( Vieira et al.ll201oh , the Atacama Cosmology Tele- 
scope at 1-2 mm, and with SPIRE (250-500 ^m; 
Griffin et all l2010n on the H e rschel Space Ob- 
servatory (IPilbratt et all l2010l: lEales et all 120101: 
Oliver et al. 2010l ) are uncovering a large number 
of these extreme SMGs, and are thus providing 
a sizable number of very bright targets that are 
ideal for follow-up CO measurements. Five of 
these extremely bright SMGs that were detected 



in the Herschel-ATLAS Science Demonstration 
Phase data were recently targeted for blind CO 
measurements with Z-Spec and Zpectrometer on 
the Green Bank Telescope (GBT), and all five 



redsh i fts were successfull y measured (|Lupu et al 
20101 iFraver et all l201ll) . These redshifts have 
been key to demonstrating that these extremely 
bright SMGs are str ongly lensed by inte rvening 
foreground galaxies (|Negrello et al 2010l ). With 
a growing number of strongly lensed SMGs being 
uncovered with Herschel/SPIRE and the SPT, 
the number of blind CO redshift measurements 
is expected to grow considerably, even before 
full operations of the Atacama Large Millime- 
ter/submillimeter Array (ALMA) are underway. 

The CO line fluxes constrain the physical prop- 
erties of molecular gas in a galaxy, including the 
total molecular gas mass as well as temperature 
and density. Since estimates of physical conditions 
depend on the relative line strengths, it is neces- 
sary to sample the full spectral line energy distri- 
bution (SLED) from the low- to high- J transitions 
in order to place useful constraints on the gas 
properties. While most galaxies at both low and 

high redshifts have only been detected in 1 — 2 CO 

1 u 11 1 1 

lines (Solomon & Vandcn Bout 2005; Grcvc et al 



20051 iTacconi et al.ll2006l l2008t iHarris et al.ll2010l: 



Ivison et al.ll2010bl : lAravena et al.ll20ldlWang etal 



2010( 1 . a growing number of galaxies have been ob- 
served in > 3 transitions in recent years. These 



both low and high redshifts ( 


Weifi et all l2005al 


Riechers et al. 2006: Ao et al. 


20081 Weifi et al. 


2007a|H; 1 


3radford et al.l 2009t Papadopoulos et al. 


2010aM_ 


Carilli et al.l l2010t iRiechers et all 201fi 


Lestrade et all 120101 iDanielson et all 201 lh as 



well as recent Herschel S PIRE /Fourier trans- 
form spectrometer (FT S; IPanuzzo et al. 2010: 



van _der Werf et alll2010h and HIFI rtLoenen et al.. 
2010h observations of nearby starbursts and active 



galaxies. For many, the full CO SLED is well de- 
scribed by a single component with warm, dense 
gas typical of star-forming regions, with kinetic 
temperatures of T^ in ~ 30 — 100 K and densities 
of riH 2 ~ 10 35 ~ 6 cm~ 3 . However, those with the 
best sampled SLEDs, spanning the full run of the 
rotational transition ladder, often require multi- 
ple gas phases to explain the observed line fluxes, 
with a warm/dense phase required to excite the 
mid- to high- J transitions, and extended, cold low- 
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excitation gas that contributes significantly to the 



J UD 5 2 line fluxes (IWard et alj|2003t ICarilli et al 



2010; iDanielson et all 120111 : IPanuzzo et al.l I20T0F 



van der Werf et all l2010t ). Several studies have 



found that most SMGs show an excess in the CO 
J = 1 —> l ine luminosity relative to the higher- J 
transitions ((Harris et al]l20ldlpraver et al.ll201lh . 
and high-resolution mapping of CO in a num- 
ber of these high-redshift galaxies reveals that 
the J = 1 — > line traces a more extended gas 
reservoir than that traced by J up > 3 transitions 
(jlvison et al.ll2010ah . 

In this paper, we present the CO redshift 
measurement and excitation modeling of HER- 
MES J105751. 1+573027 (hereafter HLSW-01), 
a multiply-imaged SMC discovered in Science 
Demonstration Phase Herschel/ 'SPIRE data as 
part of the Hersch el Multi- Tiered Ex tragalactic 
Survey (HerMES^; lOliver et al.l l2010t 2011 (in 
prep.)). High resolution imaging at 880 /im with 
the Submillimctcr Array (SMA) reveals at least 
four images separated by ~ 9", and higher res- 
olution Ap-band imaging shows strong lensing 
arcs. The lensing model (Gavazzi et al. 2011; 
hereafter Gil) suggests that HLSW-01 is lensed 
by a massive group of galaxies, with a total 
gravitational magnification of \i = 10.9 ± 0.7. 
When corrected for the magnification, this galaxy 
is found to be a very bright ultra-luminous in- 
frared galaxy (ULIRG) with an FIR luminosity of 
Lfir = (1.43 ± 0.09) x 10 13 L© and an implied 
SFR of 2460 ± lGOMgyr" 1 (Conley et al. 2011; 
hereafter Cll). The dust continuum spectral en- 
ergy distribution (SED) is warm with a dust tem- 
perature of 7d ~ 88 K determined from a simple 
one-component modified blackbody model (Cll). 
Combined with a radio flux density in excess of 
that expected from the FIR to radio correlation, 
there is tentative evidence that this source har- 
bors a bolometrically important AGN, in which 
case the SFR inferred from Lfir is overestimated. 

This paper is organized as follows: we discuss 
the Z-Spec observations and data reduction for 
HLSW-01 in Section [2] in Section [3j we describe 
the redshift determination and line flux measure- 
ments from the Z-Spec data; we describe the CO 
excitation modeling in Section and discuss the 
molecular gas properties in Section [5] We summa- 



rize these results in Section [6] We assume a flat, 
ACDM cosmology with J7m = 0.3, = 0.7, and 
Ho = 70kms~ 1 Mpc^ 1 throughout this paper. 

2. Z-Spec Observations and Data Reduc- 
tion 

We carried out the Z-Spec observations of 
HLSW-01 at the Caltech Submillimeter Observa- 
tory (CSO) from 2010 March 09 to May 12 under 
generally good to excellent observing conditions, 
with a zenith opacity at 225 GHz (monitored by 
the CSO tau meter) ranging over r 2 25 = 0.03 to 
0.1, with r 2 25 < 0.06 for 80% of the observations. 
The beamsize ranges from 25" to 40" (full width 
at half maximum) over the Z-Spec bandpass. The 
data were taken using the standard "chop-and- 
nod" mode in order to estimate and subtract the 
atmospheric signal from the raw data. The sec- 
ondary mirror was chopped on- and off-source at 
a rate of 1.6 Hz, with a chop throw of 90" while 
stepping through a 4-position nod cycle, integrat- 
ing for 20 sec at each nod position. We checked 
the pointing every 2 — 4 hrs by observing quasars 
and other bright targets located close in elevation 
to HLSW-01, making small (typically < 10") ad- 
justments to the telescope pointing model in real 
time. The total integration time (including the 
time spent in the off-source position during the 
nod cycle, but excluding all other overheads) was 
22.9 hrs. 

We analyze the data using cust omized software 



i n the same manner as described in lBradford et al 
(2009). ?oi each channel, the nods are calibrated 
and averaged together, weighted by the inverse 
square of the noise. Absolute calibration is de- 
termined by f requent (~ 1 per night) observa- 



tions of Mars (|Wrightl [20071 ) and Neptune^. The 



flux density of Mars (proportional to v 2 from 
190 — 310 GHz) decreased with apparent size from 
800 Jy to 270 Jy at 240 GHz over the course of the 
Z-Spec Spring 2010 observing season, while the 
240 GHz flux density of Neptune was 16 Jy dur- 
ing this run. We use a total of 42 planet ob- 
servations taken throughout the observing run to 
build a model of the flux conversion factor (from 
voltage to flux density) as a function of operat- 
ing ("DC") voltage for each detector separately 



1 hermes . Sussex, ac.uk 



2 http: / / smal.sma. hawaii.edu / planet vis . html 
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(|Bradford et alj|2009l) . Since the DC voltage de- 
pends on a combination of the bath temperature 
and the total optical loading on the detectors, 
we use these curves to determine appropriate cal- 
ibration factors to apply to each nod individu- 
ally. Based on the root mean square (rms) devia- 
tions of the planet measurements from the best-ht 
curves, the channel calibration uncertainties are 
4 — 9%, excluding the lowest frequencies for which 
a good model of the atmosphere is hindered by the 
183 GHz atmospheric water line. These errors are 
propagated through the data reduction and inte- 
grate down as the square root of the number of 
nods included in the average. 

Since blind CO line detections with Z-Spec re- 
quire small channel-to-channel variations, we de- 
termine small bandpass corrections to the cali- 
bration for each channel from continuum mea- 
surements of bright mm-sources. We carried out 
a total of 67 observations of bright continuum 
sources (1mm flux density, S±.i > 1 Jy), includ- 
ing J1055+018, J1637+574, J0854+201, 3C273, 
3C279, and 3C345. Each was reduced using the 
calibration determined from the planet observa- 
tions as described in the previous paragraph. We 
then fit a power-law to the continuum, separately 
for each observation, and for each channel com- 
pute a multiplicative correction factor to apply to 
the observed spectrum. We exclude the measure- 
ments from the 12 lowest frequency channels (y < 
190 GHz, contaminated by the atmospheric water 
line) and from the 232.3 GHz channel (known to 
be unstable) in the continuum fits. We then av- 
erage over all observations to determine a single 
correction factor for each channel. These bandpass 
corrections are small (ranging from 0.96 — 1.06), 
with little variation (< 5%) over all observations. 
We apply these small correction factors to the cali- 
brated data to improve the channel-to-channel cal- 
ibration for our data. 

The rms uncertainties on the final co-added 
spectrum of HLSW-01 range from 1.8 to 4.9 mJy, 
or 15% for most detectors. These errors do not 
include the uncertainties on the brightness tem- 
peratures of Mars and Neptune, which are ~ 5%. 



3. Redshift, Continuum, and Line Flux 
Measurements 

The 190-310 GHz spectrum for HLSW-01 mea- 
sured by Z-Spec is shown as a histogram in Fig- 
ure[T] The error bars represent the la photometric 
errors on the measurements (not including calibra- 
tion errors). 

3.1. Redshift Determination 

We use a custom-developed algorithm to deter- 
mine the redshift of HLSW-01 based on the de- 
tection of multiple emiss ion lines. This algo rithm 
is described in detail in ILupu et~aL (I2OIOI here- 
after L10), and we summarize it here. Since Z- 
Spec (with channel widths ranging from 720 to 
1290 km s -1 ) does not spectrally resolve the line 
emission from typical galaxies, the signal from a 
given line will approximately fall within a single 
Z-Spec channel, although the spectral response 
profiles of adjace nt channels overlap somewhat 
( Earle et al. 2006). Using a reference line list con- 
taining all CO transitions up to J = 17 — > 16 and 
the fine structure lines from neutral carbon ( [C I] ) 
and singly ionized carbon and nitrogen ([C II] and 
[N II], respectively), we compute two different es- 
timators based on combinations of the continuum- 
subtracted signal-to-noise {S/N) measured in all 
channels where the reference lines would lie for 
a given redshift. For HLSW-01, both estimators, 
Ei(z) and E%(z) (see L10 for details), are max- 
imal at z = 2.958 ± 0.007, and with maximum 
values of E 1 (2.958) = 9.1 and E 2 (2.958) = 5.3, 
the redshift is determined with » 99.99% con- 
fidence (L10). Subsequent measurements of the 
CO J = 5 — > 4 line using the Institut de Radioas- 
tronomie Millimetrique (IRAM) Plateau de Burc 
Interferometer (PdBI), the CO J = 3 -)• 2 line us- 
ing the Combined Array for Research in Millime- 
ter Astronomy (CARMA), and the CO J = 1 — s- 
line using the GBT Zpectromcter, confirmed this 
redshift (z = 2.9574±0.0001; Riechers et al. 2011; 
hereafter Rll). 

3.2. Line and Continuum Fitting 

We fit the spectrum of HLSW-01 to a model 
consisting of a power-law continuum and CO line 
emission. Since the spectrometer outputs only a 
single value per channel, it is not critically sam- 
pled, and so we use the spectral response profile 
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Fig. 1. — The mm- wavelength spectrum for HLSW-01 measured by Z-Spec (histogram). The error bars 
show the la photometric errors on the measurements and do not include the 5% uncertainty on the absolute 
flux calibration. The solid curve shows the best-fit model to the dust continuum and the CO and [C I] line 
emission, with the positions of the lines marked. This model includes the convolution of the intrinsic line 
profile (Av = 350 km s -1 ) with the spectral response profiles of the channels. 



for each channel in the fitting. We exclude in the 
fit all channels with v < 190 GHz due to poor cal- 
ibration. We fix all line widths to 350 km s -1 (full 
width at half maximum) , which is the line width 
measured by the PdBI, CARMA, and Zpectrom- 
eter assuming a Gaussian profile (Rll); however, 
we note that the fitted line fluxes are fairly in- 
sensitive to the choice of line width. The CO 
,7 = 7—5-6 line is separated from the [C I] 
3 P 2 — > 3 .Pi (hereafter [C I] 2 — > 1) fine struc- 
ture line by ~ 1000 km s -1 , or roughly one Z-Spec 
channel. For this reason we fix the redshift to 
z — 2.9574, the value measured by the PdBI, and 
include the [C I] 2 — > 1 line in the fit. The line 
fluxes (or 5cr upper limits) for the four CO lines 
in the Z-Spec bandpass and the [C I] 2 — > 1 line 
are shown in Table [TJ along with the Zpectrome- 
ter CO J = 1 -> 0, CARMA CO J = 3 ->• 2, and 
PdBI CO J = 5 -> 4 measurements (Rll). We 
were particularly unlucky in that the J = 8 — > 7 
line falls within a noisy channel of the Z-Spec 
bandpass, so we obtain only an upper limit for 
this line. The best-fit to the dust continuum is 
F v = (15.3 ±0.2) ( M ^ H -) 3 ' 2±01 m Jy. Given a 



dust temperature of — 88 K determined from 
fitting the FIR to mm continuum (Cll), we are 
not strictly in the Rayleigh- Jeans limit over the 
entire Z-Spec bandpass at z = 2.9574; however, 
this power-law provides a good fit to these data 
(reduced \ 2 01 0.97) and is sufficient for the pur- 
poses of baseline fitting. This best-fit model is 
over-plotted in Figure [TJ 

The formal errors on the line fluxes derived 
from the above fit do not include the uncertainty 
of the Z-Spec frequency scale, which is a v ~ 
lOOkrns -1 . This has important implications for 
the CO J = 7 -> 6 and [C I] 2 -> 1 line fluxes, 
as these lines are blended. We include this addi- 
tional uncertainty on the lines fluxes by shifting 
the redshift by ±a z = ±(1 + z) a v /c = ±0.001, re- 
fitting the line fluxes, and taking the upper and 
lower bounds from the statistical la errors for 
these fits. These ±lcr uncertainties in the best- 
fit fluxes are listed in parentheses in Table [TJ For 
the CO J = 7 -> 6 and [C I] 2 -> 1 lines, including 
the frequency scale uncertainty roughly doubles 
the error bounds on the measured lines fluxes. 
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Table 1 

CO and [C I] Transitions Observed in HLSW-01 



Transition 


frest 


fobs 


Flux Dcnsity a 


Line Luminosity a,b 


Line Flux a ' b 


Reference 




(GHz) 


(GHz) 


( Jy km s - 1 ) 


(10 10 K km s _1 pc 2 ) 


(1Q- 20 W m~ 2 ) 





CO ./ = 


1 - 


> 


115.27 


29.13 


1.1 ± 0.1 


4.04 ± 0.39 


0.010 ± 0.001 


GBT Zpcctromctcr 


CO J = 


3 - 


> 2 


345.80 


87.38 


9.7 ± 0.5 


3.83 ± 0.20 


0.26 ± 0.01 


CARMA; Rll 


CO J = 


5 - 


> 4 


576.27 


145.61 


23.6 ± 1.4 


3.34 ± 0.20 


1.05 ± 0.06 


IRAM PdBI; Rll 


CO J = 


7 - 


+ 6 


806.65 


203.83 


34.6 ±4.2 (±l\) 


2.50 ±0.30 (±°;™) 


2.2 ±0.3 (tit) 


Z-Spcc; this work 


CO ./ = 


8 - 


+ 7 


921.80 


232.92 


< 20 (< 30)' 


< 1.1 (< 1.7) 


< 1.5 (< 2.2) 


Z-Spcc; this work 


CO J = 


9 - 


> 8 


1036.91 


262.01 


12.7 ± 3.3 (lg 3) 


0.56 ±0.15 (±°-\l) 


1.0 ±0.3 ( +0 l) 


Z-Spcc; this work 


CO J = 


10 


-> 9 


1151.99 


291.09 


14.7 ±3.7 


0.52 ±0.13 (t° il) 


1.3 ±0.3 (t° :t) 


Z-Spcc; this work 


[C I] 3 P 


2 — 


• 3 Pi 


809.34 


204.51 


< 25 (< 48) 


< 1.8 (< 3.5) 


< 1.6 (< 3.0) 


Z-Spcc; this work 



a For the Z-Spcc data, the initial errors and 5cr upper limits come from the formal statistical errors from fitting the spectrum. The 
errors and 5cr upper limits in parentheses include the uncertainty in the Z-Spcc frequency scale as described in Scction l3.2l 

b Corrected for magnification assuming fi — 10.9. 



The CO J = 10 — >• 9 line appears some- 
what broader than the other transitions and is 
brighter than the J = 9 — >• 8 line; we there- 
fore consider possible blending with the 0-H2O 
3 12 221 water line at v Ies t = 1153.13GHz. 
Recent observations of submm-bright AGN-host 
galaxies show evidence for water vapo r emission 
([Bradford et alJl2009t ILupu et alJl2010h . the most 
striking example being the local ULIRG Mrk 
231, where the Herschel SPIRE/FTS spectrum 
reveals seven rotational emis sion lines of water 
( Gonzalez- Alfonso et al.l l2010l ) . Given the close 
proximity in frequency of the CO J = 10 — > 9 
and o-H 2 3i2 — > 2 2 i lines, we cannot deblend 
these lines within the Z-Spec spectrum; however, 
theoretical arguments supported by observed wa- 
ter line ratios in nearby galaxies suggest that the 
CO J = 10 -5- 9 line flux for HLSW-01 is not 
contaminated by water emission. The 0-H2O 
3i2 — > 2 2 i line is part of a de-excitation cas- 
cade process, following the excitation of the 32i 
level through the o-H 2 2i 2 — > 3 2 i 75 /zm transi- 
tion. The subsequent cascade results in emission 
of the o-H 2 3 2 i -> 3i2 Orcst = 1162.91GHz), 
o-H 2 3i2 -> 3 03 (^rest = 1097.37 GHz), and 
o-H 2 Q 3i2 — > 2 2 i lines (e.g ., see Figure 2 in 
Gonzalez- Alfonso et al.l 201dl ). Photon number 
conservation implies that the sum of the pho- 
tons emitted in the 0-H2O 3i2 — > 221 and 0-H2O 
3i2 — > 3o3 lines equals the number of photons 
emitted in the o-H 2 32i — > 3i2 line. Since neither 
the o-H 2 3 2 i — > 3i 2 nor the o-H 2 3i 2 — > 3o3 
lines are detected in the spectrum of HLSW-01 



at the la level, we argue that the strength of the 
0-H2O 3i2 — > 221 line is also consistent with the 
noise in our measurements. Furthermore, the ob- 
served o-H 2 3 2 i — > 3i2 and p-H 2 2 02 — > In 
(^rest = 987.92 GHz) transitions are stronger than 



the o-H 2 Q 3i2 



^21 



line for both Mrk 231 



([Gonzalez- Alfonso et al.l 120101 ) and Arp 220 (N 
Rangwala, private communication). This adds 
observational evidence against significant blend- 
ing of the CO J = 10 -> 9 line with the o-H 2 
3i 2 — > 2 2 i line in HLSW-01, since we do not detect 
the stronger water lines in the spectrum. 

The CO SLED in u nits of Wrn" 2 and in line 
lumin osity units (L' co , ISolomon fc Vanden Bout 
2005) is shown in Figure^ These values have been 
corrected for a magnification of /i = 10.9, as deter- 
mined from the lensing model of Gil. The SLED 
turns over at 6 < J 5, 8, as is t ypical for SMGs 
and q uasars ([Weifi et al.ll2007al lbl:' iDanielson et al 



20111 ); however, it is not possible to identify the 
turnover precisely given that the CO J = 7 — >• 6 
line flux may be overestimated due to blending 
with the [C I] 2 — > 1 emission line, and we have 
only an upper limit on the CO J = 8 — >• 7 line 
flux. 

4. Excitation and Radiative Transfer Mod- 
eling 

The SLED for HLSW-01 is well sampled from 
low- to high- J, and thus allows for a rigor- 
ous analysis of the CO line excitation. We use 
RADEX, a non-LTE radiative transfer code which 



7 



3.0 
2.5 
2.0 
1 .5 
1 .0 

0.5 
0.0 



£ 10 'o 



1 -Component 
2-Component (a) 
2-Component (b) 




uses an iterative escap e probability formalism 



(|van der Tak et all l2007h . 



CO J upper 

Fig. 2. — The CO spectral line energy distribution 
for HLSW-01, in line flux (top) and line luminosity 
units (bottom). The circles and upper limit (5<r) 
are from the Z-Spec measurements (this work), 
where the la error bars include the uncertainty 
from the Z-Spec frequency scale. The squares are 
from the GBT (CO J = 1 -> 0), CARMA (CO 
J = 3 -> 2), and PdBI (CO J = 5 -> 4) obser- 
vations (Rll). All data have been corrected for 
magnification assuming = 10.9. The black solid 
curve represents the maximum likelihood single- 
component model discussed in Section 15.11 The 
red dot-dashed curve shows the maximum likeli- 
hood solution from a two-component gas model 
described in Section 15.21 The red dashed and 
dotted curves show the contribution from the 
warm and cold components, respectively. The 
blue dot-dashed curve demonstrates a different 
two-component gas model (also described in Sec- 
tion I5.2[) that provides a good fit to these data; 
in this case, both components (blue dashed and 
dotted curves) are relatively warm. 



to model the CO ex- 
For a specified gas den- 
and CO 



citation of HLSW-01 
sity (tt,h 2 ), kinetic temperature (Tki n 
column density per unit line width (Nco/Av), 
RADEX calculates the excitation temperatures, 
line optical depths, and line surface brightnesses. 
For our analysis, we use the escape probability 
in an expanding spherical cloud and assume a 
Tcmb(^) = 2.73 (1 + z)K blackbody for the back- 
ground radiation field. However, we note that the 
RADEX results are not sensitive to the precise 
form of escape probability chosen. 

Following a sim ilar analys is to that described in 
Ward et alj (l2003l) (see also iBradford et all [2009; 



Navlor et al.ll2010l and Kamenetzky et al. 2011) 
we use RADEX to compute the expected CO line 
intensities over parameter space in T^ n , nn 2 , and 
A^cOj assuming a line width of Av = 350 km s -1 
(Rll) for all transitions. Rather than computing 
line intensities over a grid in parameter space, 
we use a Markov-cha i n Monte Carl o method 
((Metropolis et al.lll953t lHastingslll970l ) to deter- 
mine the likelihood distributions for these param- 
eters by comparing the RADEX results to the ob- 
served line intensities. The integrated line fluxes, 
Sco^v (in Jykms^ 1 , Table [T]) are converted into 
Ray leigh- Jeans equivalent velocity-integrated in- 
tensities, TrAu (in Kkms -1 ) for comparison with 
the output from RADEX: 



TrAv 



5 co Aw 



'obs 



(1 + z), 



(1) 



where /i = 10.9 (Gil) is the gravitational lens- 
ing magnification factor, and f2 cm is the intrinsic 
solid angle of the CO emitting region (see the next 
paragraph). We propagate the line flux errors us- 
ing the ltr error bounds that include the uncer- 
tainty in the Z-Spec frequency scale. We do not 
include errors in the absolute flux calibration as 
these are small compared to the random errors on 
the measurements; nor do we account for system- 
atic calibration errors among the measurements, 
as these are not possible to characterize without 
at least two measurements of the same CO tran- 
sition using different telescopes. 

The CO emission is only marginally resolved in 
both the IRAM PdBI and CARMA maps. Given 
the uncertainties in the size of the CO emitting 
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region and the gravitational magnification, com- 
bined with the fact that the area filling factor of 
the molecular gas may be < 1, we choose to in- 
clude the intrinsic source solid angle, O em , as a 
fourth parameter in our likelihood analysis. The 
intrinsic solid angle refers to the observed solid an- 
gle in the absence of lensing. Finally, we impose an 
upper limit for the kinetic temperature of 3000 K; 
above this temperature, collisional dissociation of 
CO becomes important (with a weak dependence 
on the gas density). 

We ran a 150, 000-step Markov chain to deter- 
mine the likelihood distributions of the gas pa- 
rameters for HLSW-01. The marginalized distri- 
butions for the four parameters are shown in Fig- 
ure [3] (gray, unshaded histogram) , with the maxi- 
mum likelihood values marked with gray solid ver- 
tical lines. Note that since these are marginalized 
distributions, the peaks do not necessarily coin- 
cide with the 4-dimensional maximum likelihood 
solution. We use these 1-dimensional distributions 
only for computing the marginalized 68.3% confi- 
dence interval on each of the parameters, which 
is indicated by the gray vertical dotted lines in 
Figure |3J These results are summarized in Ta- 
ble [5J Figure S] shows the 2-dimensional marginal- 
ized distributions for Tki n and 71h 2 , £l C m and Nqo , 
and Nqo and tih 2 (gray contours). For T^ n and 
njj 2 , the contours fall along a line of constant pres- 
sure. 

We also derive the gas pressure (P = uh 2 Tjd n ) 
and mass, where the latter is estimated as: 



AT™ = 



» em DlNco(lAm Sa ) 
Xco 



(2) 



where fl om D\ is the area of the emitting region, 
Da is the angular diameter distance, Xqo = 
nco/nn 2 is the relative abundance of CO to H2, 
and mn 2 is the mass of the hydrogen molecule 
(with a factor of 1.4 to account for helium). We as- 
sume Xco — 2 x 10~ 4 , a reasonable value based on 
observations of nearby ga laxies and star-fo rming 
regions m the Galaxy fe.g. lWang et~aLll2004 . The 
marginalized distributions for these secondary pa- 
rameters are also shown in Figure El 

The maximum likelihood solution gives Tk; n ~ 
570K and nn 2 ~ 3 x 10 2 cm~ 3 , a condition that 
is unlike the average molecular gas seen in most 
galaxies. This solution is also inconsistent with a 



gas that has at least enough velocity dispersion to 
correspond to virialized moti on under its own self 
gravit y. Using the notation of lPapadopoulos et al 



(2007) and reforming their Equation 2 for our pa- 
rameterization yields the following constraint for 
a self-gravitating gas: 



K vir = 



_ (dv/dr) obs 



Av nj/ 2 Xco 



(dv/dr) vil (47rGa(1.4 m H2 )/3) V* N CO 

(3) 

where a = 1 — 2.5 depending on the cloud den- 
sity profile, and we estimate the velocity gra- 
dient for an individual cloud as {dv / dr) Q \, s = 
( ■£»-) n H2 ^co (Kamenetzky et al. 2011). As- 



> 1 



\7fc~o, 

suming Xco = 2 x 10" 4 , K vil = 0.2 - 0.4 for 
the maximum likelihood solution. For this rea- 
son, we add a prior to our likelihood analysis 
that constrains parameter space to solutions where 
the gas is self-gravitating, i.e. where Equation [3] 
holds, and we repeat our Markov chain calcula- 
tions. These results are presented in Tabled Fig- 
ure [3] (shaded histograms) , and Figure 0] (black 
contours), where the maximum likelihood solution 
now gives T kin « 230K and n H2 ~ 1.2 x 10 3 cm~ 3 . 
Throughout the rest of this paper, we discuss only 
the solution where the constraint on K v { r has been 
enforced. 

5. Results and Discussion 

Our likelihood analysis for HLSW-01 sug- 
gests that the CO gas is primarily tracing warm, 
moderate-density gas with Tki n = 227 



+8 K 
141 1V 



and 



«H 2 



_ 1 9+2.3 



x 10 3 cm 



This single- 



component model provides a good fit to our 
data, agreeing within 1.5er of all measurements 
as demonstrated in Figure [3J where the model 
line fluxes corresponding to our maximum like- 
lihood solution (solid black curve) are com- 
pared to our data. However, this model falls 
short of the J = 10 — > 9 line flux, and it is 
likely that the full CO SLED is really a com- 
posite of multiple gas components, as seen in 



dWard et al. 


l2003t IWeifi et al. 


l2007al: iGreve et al. 


200a Panuzzo etaTJ l2010t lLoenen et al. 


20101 


Riechers et al.l2O10L Danielson et al.l2011: Carilli et al. 


2010: van der Werf et all 201(*. In Section IB31 



we discuss the implications our single-component 
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Fig. 3. — Top left, top center, bottom left, bottom center: Marginalized distributions for the four parameters 
from the excitation modeling of HLSW-01. Top right, bottom right: Marginalized distributions for pressure 
and gas mass, derived from the four primary parameters. All of the distributions are normalized at their 
peaks. In each panel, the gray, unshaded histogram shows the resulting distribution when the gas is not 
required to be self-gravitating, while the shaded histogram shows the distribution when this constraint is 
enforced (i.e. Equation [3]) . The vertical solid and dotted lines indicate the maximum likelihood solutions 
and the marginalized 68.3% confidence intervals, respectively, with (black) and without (gray) the constraint 
of self-gravitating gas. The source solid angle and gas mass distributions assume a gravitational lensing 
magnification factor of [i = 10.9 for HLSW-01. For a different value of [i, multiply the x-axes in the bottom 
center and right plots by (^p)- 
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Table 2 

Results of Single-phase CO Excitation Modeling for HLSW-01 



K V { T unconstrained K v ir > 1 

Parameter Maximum Likelihood Solution Maximum Likelihood Solution Units 



T kin 
riH 2 
Nco 



566 (246 - 845 ) 
0.3 (0.2 - 0.7 ) 



2.7 (1 
5.0 (4 



2.9 
5.6 



227 (86 - 235 
1.2 (1.1 - 3.5 
2.1 (1.6 - 2.5 
4.8 (4.6 - 5.9 



10 3 cm" 1 
10 20 cm" 2 

(^) io- 13 



P 

M gaB 
dv/dr 



1.6 (1 
3.9 (2 
0.22 (0.18 



5 - 2.5 ) 
i - 4.2 ) 



0.55 



2.6 (2.4 
2.9 (2.3 
1.2 (1.1 - 3.0 



3.7 ) 
3.9 ) 



10 a K cm" 

(^^H 1 ^) 1010 M ° 

km s" 1 pc" 1 



"CO 



Note. — The columns are as follows: 1) the parameter; 2) the maximum likelihood value of the parameter and 
68.3% confidence interval (in parenthesis), without the constraint that the gas is self-gravitating; 3) the maximum 
likelihood value of the parameter and 68.3% confidence interval when the gas is required to be self gravitating 
(Equation[3j; and 4) the units for the given values. 




Fig. 4. — The 2-dimensional marginalized distributions for Tkin and 7Jh 2 (left), Q cm and iVco (center), 
and Nqo and «,h 2 (right). In all plots, the contours are 68.3%, 95.5%, and 99.7% assuming Gaussian 
distributed errors. The gray contours show the resulting distributions when the gas is not constrained to 
be self-gravitating, and the gray plus symbols mark the maximum likelihood solution in that case. The 
black contours show the distributions with the constraint that the gas is self-gravitating, with the black 
stars marking the maximum likelihood solution. In the left plot, the dotted lines indicate lines of constant 
pressure (in units of K cm" 3 ), and the dashed pressure line passes through the maximum likelihood solution. 
In the center plot, the dotted lines show constant gas mass (in units of M ), and the dashed line passes 
through the maximum likelihood solution. The dotted lines in the right plot show constant velocity gradient 
(in units of kms" 1 pc" 1 ), where the dashed line indicates the maximum likelihood solution. The source solid 
angle assumes a gravitational lensing magnification of fi = 10.9 for HLSW-01. For a different value of //, 
multiply the j/-axes in the center plot by ( i2Ji J . 
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gas model would have for HLSW-01, and in Sec- 
tion 15.21 we consider a few two-phase gas models 
that are also consistent with our data. 

5.1. Single-component Gas Model 

For our best-fit single-component gas model, all 
of the CO transitions are sub-thermally populated, 
with an excitation temperature of T ex ~ 135 K for 
the CO J = 1 — > line, and decreasing to T cx ~ 
25 K for the CO J = 10 -4 9 line. The high gas 
kinetic temperature (1km = 230 K) is required to 
excite the higher- J transitions, and is considerably 
larger than the dust temperature of Td = 88 K 
determined from Cll. This suggests that heating 
mechanisms that effectively transfer energy into 
the gas (e.g. X-ray dominated regions, or heating 
by the dissipation of turbulence or cosmic rays) 
may be present in this system (although colder 
gas solutions that match the dust temperature are 
also consistent with our data). 

Our result for HLSW-01 is similar to that 
for the local starburst galaxy M82, where the 
■Aip > 4 lines are best-fit by warm gas with mod- 

i 5 x 10 3 cm~ 3 , 
in contrast to 
M82, this warm gas is able to account for the 
excitation of the lower- J lines as well, whereas in 
M82 the J up < 3 line fluxes are dominated by 



erate density (Jkin ~ 550 K, nn 2 ; 

" 2010h . However, 



Panuzzo et al 



Msz tne J up S o line nuxes are dominate d by 
a cold, less dense component Jwfild et al.l 1992; 



Guesten et al.lll993t [Mao et al.ll2000t IWard et al 



20031: IWeifi et al.l l2005bf ). Given its low critical 
density, the CO J = 1 — > line is easily excited 
in both cold, diffuse gas as well as warm, dense 
gas in star-forming regions. There is growing ev- 
idence that many high-redshift SMGs contain a 
significant amount (^50% of the total mass) of 
extended, cold gas that is not associated with 



the nuclear starburst (e. 



Ivison et al.. 
2010) 



2010a: 



Danielson et al.1 120111: iCarilli et al.l 12010( 1 . This 
leads to an excess in the CO J = 1 — > line 
luminosity relative to the higher- J transitions, 
with a line ratio of i?3 ; i = L ' C o ( ^_ 9 ^/ L' ran _ r[ ^ ~ 
0.6 le.g. lHarris et all l2010t iFrayer et alJ 12011 



Ivison et al.l l2010al l and a line wi dth ratio of 
Af ca n -o) / A«co(3-2) = 1.15 ±0.06 (jlvison et al. 
2010al ). However, this does not appear to be the 
case for HLSW-01, where R3.1 — 0.95 and the ob- 
served line widths for the J up = 1, 3, and 5 transi- 
tions are all consistent with the presence of a single 
component. This suggests that the J up < 5 tran- 



sitions for HLSW-01 are largely tracing the same 
volume of molecular gas in this galaxy, and given 
the large uncertainties in the higher- J observa- 
tions, we see no strong evidence that a two-phase 
gas model is needed to explain the CO excitation 
for this source. 



We derive 
4_a+ 



-1.1 / iojA 
'0.2 \ n ) 



an equivalent diameter of 1.25 +0 ! 



a source solid angle of Q em = 
x 10~ 13 sr, which corresponds to 

0.03 ) k P C 

at z = 2.9574, assuming spherical geometry. 
Combined with the warm temperature and mod- 
est density, this suggests that the molecular gas 
is excited by an intense starburst that is less 
centrally concentrated than that seen in local 
ULIRGs, but consistent with previous findings of 
exte nded star formation i n SMGs (1.5 — 3 kpc , 



e.g. iBiggs fc 



Lestrade et al 



vison 



20081 ITacconi et al.l I2008L 



2010t ISwinbank et al.1 120T0I1 . It 



is possible that the molecular gas re servoir is ex- 
tended over an even large r area (e.g.. Ilvison et al 
2010al : ICarilli et alj|2010h . and tt cm really repre- 



sents the product of this area with a filling frac- 
tion. The velocity gradients observed for the four 
images in the CO J = 3 ->• 2 and CO J = 5 -)• 4 
maps suggest a complex velocity structure, possi- 
bly arising from a major merger (Rll); however, 
higher spatial resolution CO maps are needed in 
order to confirm this interpretation. 

The total molecular gas mass derived from 
this line excitation modeling (Equation [2]) is 

2 - 9 ±ae (^) ( 1 7r) >< ^°M Q , similar to 
those of other high - redshi ft SMGs and quasars 
fe.g. ITacconi et al.l l2008t ICoppin et al.l l2008t 
' " 20101: iRiechers et al.l |2010| ). This 



Wang et al 



value is consistent with the estimate from ^co(i-o) 
(3.3 x 1O 1O M0, Rll) assuming a conversion fac- 
tor of aco = 0.8 Mq (K km s _1 pc 2 ) -1 , a value 
appropriate for ultra-luminous infrared galaxies 



(Downes fc Solomonll998l:ISolomon fc Vanden Boud 



20051 ). Cll derive an SFR of 2460 ± 160 A/q yr" 1 
(corrected for le nsing) from the FIR luminosity 
( Kennicutt 19981) . Assuming that all of the molec- 
ular gas will be consumed by star formation, this 
implies a gas depletion times cale of lUj^' 2 , Myr 



a value similar to ( 12 Mvr. Ilvison et al. ||2010a) 



or somewhat shorter than (~ 40 Myr 



2005; ISolomon fc Vanden Boutl 120051 1 the mean 



Greve et al 



gas depletion timescale of typical SMGs with sim- 
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ilar luminosities. However, it is possible that the 
SFR is overestimated if a significant fraction of the 
dust is heated by AGN activity, and this estimate 
also neglects feedback mechanisms. 

At the resolution of the existing interferometry 
measurements, the SMA 880 /im and CO J = 3 — >• 
2 and J = 5 — > 4 maps show consistent source 
sizes, suggesting that the dust and gas are well 
mixed and therefore trace the same regions. We 
estimate the dust mass from the optical depth and 
the dust absorption coefficient, where 



M, 



dust 



— o n 2 — 

a dorm ' J A i 



(4) 



k v = 2.64 m 2 kg -1 0/2400GHz)^ (jDunne et al 
20Q3h . and r„ = (v/v ) p , with p = 1.94 ± 0.14 
and vq = 1550 ± 150 GHz determined from fit- 
ting a modified blackbody to the full FIR-mm dust 
SED of HLSW-01 (Cll). We derive a dust mass 



of M d = 5.2 



+ 1.6 
-1.1 



10 Mq (not including errors 



on k„), and a dust to gas mass ratio of < 1/55, 
considering that the total gas mass could be larger 
since we only estimate the mass of the molecular 
gas. This is consistent with the dus t to gas mass 
ratio s of typical SMGs (~ 1/60; e.g. Kovacs et al 



2006 : ICoppin et al.l l2008t iMichalowski et al.l l2010: 



Santini et al.l 2010h . Our dust mass estimate is 
~ 5 times higher than that estimated via £fir 
from the single-component dust model of Cll, as- 
suming the same k v . This suggests that 1) a 
single-component model for the dust (and/or CO) 
emission is not a good representation of the physi- 
cal properties of this galaxy, and/or 2) the gas and 
dust are not co-spatial. Longer- wavelength contin- 
uum measurements are needed to detect the pres- 
ence of multi-component dust, while higher spatial 
resolution imaging of the dust and CO would allow 
us to determine their relative spatial distribution. 

Recent results by IPapadopoulos et al. (|2010bl) 
sugg est that for many (U)LIR Gs, including Arp 
220 (jPapadopoulos et al 1 l2010ah . high dust extinc- 
tion can suppress the high- J CO line fluxes. As- 
suming a single-component dust model for HLSW- 
01 (Cll), the dust optical depths at the CO 
J up = 7—10 rest frequencies are t v ~ 0.3 — 0.6, 
and thus may be non-negligible. However, we can- 
not definitively know whether this effect is impor- 
tant for HLSW-01 given the degeneracy between 
high dust optical depths and large amounts of cold 
dust. A two-temperature component fit to the 



dust SED is poorly constrained and does not rep- 
resent an improvement over a single-component 
model (Cll). However, this fit only includes rest- 
frame frequencies v > 900 GHz, and it is possible 
that a significant amount of col d dust which i s best 

2000; 



traced at lower freq 



uencies ( Dunne et al 



Dunne fc Ealesl 120011 ) is present in this galaxy. 
Furthermore, to properly include this effect in the 
excitation modeling would require knowledge of 
the relative geometry of the molecular gas and 
dust, which cannot be deduced from our current 
data. High resolution mapping of the molecular 
gas and dust emission may shed light on this issue 
in the future. 

5.2. Two-component Gas Models 

The single-component gas model discussed in 
the previous section provides a good fit to the CO 
SLED for HLSW-01. However, in light of recent 
results showing that the CO excitation for a large 
number of local starbursts and high-redshift SMGs 
requires multiple gas phases, we explore a few two- 
component gas models in this section. With only 
six line fluxes measured, we cannot well constrain 
a full 8-parametcr fit; however, we can place sim- 
ple constraints on some of the parameters in order 
to test solutions that resemble those of other high- 
redshift SMGs. 

We explore the possibility that the CO ex- 
citation for HLSW-01 arises from the combi- 
nation of a dense, star-forming component and 
cold, lower-density gas that is not necessarily as- 
sociated with star-formation. We require that 
the star-forming component has a gas density of 



«H 2 



10 



3.5-6 



cm and a kinetic temperature 



of Tkm > 30 K, which is typical of star-forming 



regio n s in starburst nu c lei (e.g. iRiechers et al 



20101 ICarilli et all l2010l IPanielson et al.l l201lh 



We then require that the other component has 
a lower density and temperature than the star- 
forming component, but a larger emitting region 
(or alternatively, filling fraction). We use the 
same prior on the velocity gradient described in 
Section |4] for the star- forming component, but we 
relax the prior on the velocity gradient of the sec- 
ond component so that it is not required to be 
self-gravitating. The 2-dimensional marginalized 
distributions for T^ n and nn 2 for each component 
are shown in Figure [5] 

The maximum likelihood solution for this two- 
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Fig. 5. — The 2-dimensional marginalized distri- 
butions for Tki n and nn 2 for each component in 
the two-phase gas model described in Section 15.21 
The contours are 68.3%, 95.5%, and 99.7% as- 
suming Gaussian distributed errors, and the stars 
mark the maximum likelihood solution (model "2- 
Component (a)" in Figure [2]). The crosses mark 
another solution that provides a good fit to the 
data (model "2-Component (b)" in Figure [5]). 



component model under the above constraints is 
T kin = 180 K, tih 2 = 1.6 x 10 5 cm- 3 for the star- 
forming component, and T^ in = 30 K, nu 2 = 
8.1 x 10 4 cm~ 3 for the other component. This 
model is shown in Figure^ (model "2-Component 
(a)"), where the dashed red curve corresponds to 
the star- forming component, the dotted red curve 
corresponds to the cooler component, and the dot- 
dashed red curve shows the summed SLED from 
both components. The best-fit model thus consists 
of two dense gas phases, with warm gas distributed 
over a small (330 pc equivalent diameter) region 
(assuming a filling fraction of 1), and cold gas dis- 
tributed over a larger (2.4 kpc diameter) region. 
The total molecular gas mass is 5.1 x 1O 1O M0, 
with 30% and 70% of the gas mass in the warm and 
cold phases, respectively. The warm, dense com- 
ponent remains thermalized out to J up = 10, and 
provides the best match to our J up = 9 — 10 mea- 
surements, where these lines have roughly equal 
brightness temperatures. 

However, the likelihood space for this two- 
component model does not have a single well 
defined peak, and many other solutions repro- 
duce the lines fluxes nearly as well. For exam- 
ple, a warm/dense component with T^m = 210 K, 
n H2 = 5. 



10 cm , and a warm/moderate- 



density component with X^i, 



130 K, riH 2 



1.3 x 10 cm -3 provides a good fit to the data as 
well. This model is shown by the blue curves in 
Figure [5] (model "2-Component (b)"), where the 
dashed, dotted, and dot-dashed curves show the 
warm/dense and warm/moderate-density compo- 
nents, and their sum, respectively. In this sce- 
nario, the warm/dense gas is concentrated in 
a region ~ 270 pc, while the warm/moderate- 
density gas is extended over a larger region 
(~ 1.4 kpc). The total molecular gas mass is 
4.3 x 10 10 M Q , with 15% and 85% of the gas mass 
in the warm/dense and warm/moderate-density 
component, respectively. We note that this lat- 
ter component is similar to the results from our 
one-component fit, where the addition of a small 
fraction of warm/dense gas brings the J = 10 — > 9 
model flux in better agreement with the data. 

Our data show no evidence for two-component 
models with a significant contribution from ex- 
tended, cold, moderate-density gas with Tfa D ~ 
30 - 5 K and mr-, ~ 10 2 ~ 3 ,5 cm -3 like that seen in 
M82 ( Panuzzo et al.ll2010l ) and many high-redshift 
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galax i es (e.g Riechers et al. 2010t Danielson et al 



20111 : iCarilli et al.l I2010D . as can be seen in Fi. 
ure El While recent resu lts bv lHarris et al. ( 201 
and Ivison et al. I (l2010ah suggest that most SMGs 
contain a large fraction (> 50% of the total) of ex- 
tended, cold gas with moderate-density, HLSW-01 
appears to be an exception. 

5.3. Excitation Mechanisms 

The CO measurements for HLSW-01 do not 
have the statistical power to discriminate between 
various mechanisms for heating the molecular gas, 
such as UV radiation from star formation in pho- 
ton dominated regions (PDRs) or X-rays from an 
AGN in an X-ray dominated region (XDR). The 
Jup > 7 line fluxes measured with Z-Spec have 
large uncertainties, and it is these high-J transi- 
tions that can most strongly discriminate between 
such heating mechanisms. While high gas temper- 
atures like that found for HLSW-01 (for both the 
single- and two-component models) are usually as- 
sociated with heating in XDRs for sources that are 
known to host powerful AGNs (jWeifi et al.ll2007a : 
Ao et al.ll2008l ). heating by cosmic rays or dissipa- 
tion of turbulence can also boost the high-J line 
fluxes; these latter mechanisms are thought to be 
responsible for the high exci tation for M82, whic h 
does not host a strong AGN ( Panuzzo et al.ll201oh . 
Based on the high dust temperature and 1.4 GHz 
radio emission (in excess of that expected from 
the FIR-radio relation, Cll), there is tentative 
evidence for a bolometrically important AGN in 
HLSW-01. 

Obtaining higher S/N measurements of the 
Jup = 7 — 10 transitions - which are all accessi- 
ble from ground-based telescopes - and spectrally 
resolving these lines would go a long way to help 
differentiate between various one- and two-phase 
gas models, and to better constrain the heating 
mechanism for the gas. In particular, high spectral 
resolution is needed to properly deblend the CO 
J = 7 ->• 6 line from the [C I] 2 — > 1 line. Directly 
measuring the line widths for these high-J transi- 
tions can also help to determine whether they are 
tracing the same regions as the low- J lines. 

6. Conclusions 

We have measured a CO redshift of z = 
2.958 ± 0.007 for the multiply-imaged, lensed 



galaxy HLSW-01 using Z-Spec on the CSO. Four 
rotational transitions of CO lie within the Z-Spec 
bandpass ( J up = 7 — 10) at this redshift, and we 
have used a well tested algorithm (L10) to robustly 
identify the redshift for this source with ^> 99.99% 
confidence. This redshift has been confirmed with 
subsequent observations with the IRAM PdBI, 
CARMA, and the GBT Zpectrometer. 

With six CO line measurements and one upper 
limit, the SLED for HLSW-01 is well sampled from 
low- to high-J, and we carry out a likelihood anal- 
ysis using predictions from RADEX to constrain 
the molecular gas properties for this galaxy. We 
find that a single-component, warm/moderate- 
density gas with Tkin = 86 — 235 K and nn 2 = 
(1.1 — 3.5) x 10 3 cm~ 3 provides a good fit to these 
data. However, several two-component gas mod- 
els can describe the CO SLED as well; higher S/N 
observations of the high-J lines with finer spectral 
resolution are needed in order to differentiate be- 
tween these scenarios. 

Based on the nearly equal line luminosities mea- 
sured for the J = 1 — > 0, and J = 3 — > 2 tran- 
sitions, these data are inconsistent with models 
that require an extended, cold/moderate-density 
gas component. This is in contrast to recent re- 
sults that show that the CO J = 1 — ► line flux for 
most high-redshift SMGs, as well as nearby star- 
burst galaxies such as M82, requires the existence 
of such a low-excitation component. High spatial 
resolution CO maps of both low- and high-J lines 
with CARMA, the IRAM PdBI, and/or the EVLA 
could potentially provide more concrete informa- 
tion on the presence and distribution of cold and 
warm molecular gas within this galaxy. 
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